Tomography and spectroscopy as quantum computations 
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Determining the state of a system and measuring properties of its evolution are two of the most 
important tasks a physicist faces. For the first purpose one can use tomography, 1 a method that after 
subjecting the system to a number of experiments determines all independent elements of the density 
matrix. For the second task, one can resort to spectroscopy, a set of techniques used to determine 
the spectrum of eigenvalues of the evolution operator. In this letter, we show that tomography and 
spectroscopy can be naturally interpreted as dual forms of quantum computation. We show how to 
adapt the simplest case of the well-known phase estimation quantum algorithm to perform both tasks, 
giving it a natural interpretation as a simulated scattering experiment. We show how this algorithm 
can be used to implement an interesting form of tomography by performing a direct measurement 
of the Wigner function (a phase space distribution) of a quantum system. We present results of such 
measurements performed on a system of three qubits using liquid state NMR quantum computation 
techniques in a sample of trichloroethylene. Remarkable analogies with other experiments are dis- 
cussed. 

The basic idea discussed in this letter can be described in terms of the quantum algorithm represented 
by the circuit shown in Figure 1 (see |2| for an introduction to quantum circuits and algorithms): A system, 
initially in the state p, is brought in contact with an ancillary qubit prepared in the state |0). This ancilla 
acts as a "probe particle" in a scattering experiment. The algorithm is: i) Apply an Hadamard transform to 
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the ancillary qubit (where H\0) = (|0) + |l»/\/2, H\l) = (|0) - \l))/y/2), ii) Apply a "controlled-*/" 
operator (does nothing if the state of the ancilla is |0) but applies U to the system if the ancilla is in state 1 1)), 
iii) Apply another Hadamard gate to the ancilla and perform a weak measurement on this qubit detecting its 
polarization (i.e., measuring the expectation values of Pauli operators a z and a x ). The above circuit has the 
following remarkable property: 

(a z ) = Re[Tr(0» ], (a x ) = -Jm[Tt(Up) }. (1) 

Thus, the final polarization measurement reveals a property determined both by the initial state p and the 
unitary operator U. Versions of the above idea play an important role in many quantum algorithms. For 
pure input states, it occurs in Kitaev's solution to the Abelian stabilizer problem (a generalization of the 
factoring problem). 3 This was adapted by Cleve et a/Pto revisit most quantum algorithms, with the circuit 
of Figure 1 being the simplest instance (see [ 5 1 for yet another presentation of the algorithm as a tool for 
physics simulations). The extension to mixed states and a version of Equation 1 is in [6 1. 

Here, we observe that as a result of Equation 1, this circuit can be used for dual purposes: We can use 
it to extract information on the operator U if we know the state p. Alternatively, we can use it to learn 
about the state p by using some specific operators for U. In this sense, this circuit can be adapted either as a 
tomographer or as a spectrometer which, therefore, are dual faces of this quantum computation. Moreover, 
it is interesting to think of the above algorithm as simulating a scattering experiment: The ancillary qubit 
plays the role of a probe particle: it interacts with the scatterer (the system, initially in state p) and is later 
detected. From the statistics gathered by the detector measuring the polarization of the probe, we can either 
reconstruct the state of the scatterer (if we know about the interaction) or learn about the interaction (if we 
know about the state of the scatterer). Describing the algorithm in physics language helps in establishing 
analogies between the above proposal and some already performed experiments that, interestingly enough, 
can be interpreted as special instances of the algorithm (see below). 

Let us first discuss how to use this to build a spectrometer. The idea is simple: For example, prepar- 
ing the system in a completely mixed initial state with p = I/N (N is the dimensionality of the space 
of states of the system), the final measurement is (a z ) = Re[Tr(C7) }/N. This is directly proportional to 
the trace of U, the sum of all its eigenvalues, and therefore has information about the spectrum. More- 
over, we can build a real spectrometer by simply adding an extra register (with n\ qubits) and a controlled 
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Fourier transform. With this new circuit, shown in Figure 2, the final polarization is (a z ) = g(2E) = 
Re[ Ya exp(4vri£t/T) Tr(J7*) }/NT (where T = 2™ 1 - 1). The function g(E) is the Fourier transform of 
Tr(J7*), which is nothing but the spectral density of U (smoothed over a scale 2-k/T). One can use this same 
idea to design circuits to measure, for example, the structure function, defined as the Fourier transform of 
|TrC/*| 2 , characterizing the spectral correlations and the level spacing statistics of U, which are of interest 
in studies of quantum chaos 7 (some other related constructions and efficient networks can be found in [ 8 1). 
It is also important to notice that the above is a rather general spectrometer since it can be adapted to study 
properties of non-Unitary operators. For example, by adding an extra ancillary register one can use it as 
a spectrometer for any operator obtained as a weighted sum of unitary operators. We emphasize that the 
above is a simulation of a spectrometer where one has complete control over the operator U that acts on 
the system. This is the case for any physics simulation, either quantum or classical, where our simulation 
is intended to be used for the purpose of exploring models. On classical computers, the task of providing 
spectral information as accomplished by this quantum algorithm appears to be exponentially hard in gen- 
eral. Of course, complete determination of the spectrum of U is still inefficient, unless the eigenvalues are 
hugely degenerate. But using our simulated spectrometer one can efficiently determine important specific 
properties of the spectrum, such as directly sampling the spectral density or the structure function. 

The circuit in Figure 1 can be adapted for another important purpose: state tomography. Every time 
we run the algorithm for a known operator U, we extract information about the state p. Doing so for a 
complete basis of operators {A(a), a = 1, . . . , N 2 — 1} one gets complete information and determines the 
full density matrix. Different tomographic schemes are characterized by the basis of operators A(a) one 
uses. In NMR spectroscopy, for example, a convenient basis set is formed by all tensor products of Pauli 
operators (the basis of the product operator expansion 9 ). This kind of tomography is particularly well suited 
in this context since all A{ot) are easily reached from observables by a sequence of r.f. pulses and periods 
of free evolution.^ However, as the above is a generic tomographer one can explore other choices. In 
particular, we can use it to implement another important tomographic scheme char acterized by the operators 

A(q,p) = ^U q RV~ p exp(i2TT P q/2N), (2) 

where U is the operator producing a cyclic shift in the computational basis (U\q) = \q + 1)), V is the shift in 
the basis related to the computational one via the discrete Fourier transform and R is the reflection operator 
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(R\q) = \N — q)). By doing this, we measure the discrete Wigner function of the state of the system, 
which is defined as M^(q) = Tr(A(a)p) (we use a = (q,p)). This function is the basic tool to represent 
the state of a quantum system in phase space, the natural arena of classical physic^ and has often been 
used to understand the nature of the quantum to classical transition.^! Its use is common in various areas 
of physics, but has been mostly restricted to continuous systems (discrete Wigner functions have been used 
only recently in the context of quantum computation 1 ^). The defining properties of the Wigner function 
are: (PI) W(a) is real valued, (P2) Inner products between states pi and p% can be obtained as Tv(p\p2) = 
N Wi(a)W2(ct), (P3) Adding the value of W(a) over all points in any phase space line gives a positive 
number which is the probability of measuring a physical observable. These are consequences of interesting 
properties of the phase space point operators j4(a).U2l2l Q ne Q f th ese properties is the completeness of the 
set A(a) (a complete orthogonal set is obtained with a restricted to the first N x N sub-grid of the phase 
space but it is necessary to define W(a) in a grid of 2N x 2N points for P3 to hold 1 ?* 14 !). One can expand 
the state p in terms of this set as p = N J2 a W{a)A(a). The coefficients of the expansion are the Wigner 
function of p. 

Every time one runs this algorithm with U = A(a) a direct measurement of W(a) is obtained (for any 
given a). Before presenting results of such measurement for a simple system, it is worth mentioning that 
the same scheme can be used to measure W{a) in more general cases. For a continuous system W(a) is 
defined as the expectation value of A(a) = (a)RD{a) / irk where D(a) is the phase space displacement 
operator D(a) = e~x.p(i(pQ — qP)/h). Therefore, to measure W{a) we need to run the algorithm with the 
system in the state D(a)pD^ (a) (obtained simply by displacing p) and use U = R} 5 Measuring directly 
the Wigner function has been the goal of a series of experiments in various areas of physics (all dealing with 
continuous system*^). It is remarkable that our measurement scheme describes the recent experiment that 
determined TV (a) for the electromagnetic field in a cavity QED setup jilliU In this case the system is the 
mode of the field stored in a high-Q cavity and the ancillary qubit is a two level atom. The measurement of 
the Wigner function is done via a scattering experiment that directly follows the steps described in Figure 1 : 
i) The atom goes through a Ramsey zone that has the effect of implementing an Hadamard transform. An 
r.f. source is connected to the cavity displacing the field (by an amount parametrized by a) and preparing 
the state D(a)pD^ (a); ii) The atom goes through the cavity interacting dispersively with the field in such 
a way that only if the atom is in state | e) does it acquire a phase shift of it per each photon in the cavity 
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(i.e., this interaction is a controlled-exp(— iirN) gate, where N is the photon number, which is nothing but 
a controlled reflection), iii) The atom leaves the cavity entering a new Ramsey zone and is finally detected 
in a counter either in the \g) or |e) state. The Wigner function is measured as the difference between 
both probabilities: W(q,p) = 2(P(e) — P{g))/h. As we see, this cavity-QED experiment is a concrete 
realization of the general tomographic scheme described above. 

For a system with a finite dimensional Hilbert space the algorithm can be efficiently decomposed as 
a sequence of elementary steps: All controlled- U , V and R operations can be implemented via efficient 
networks like the ones shown in 1 19 1 which require a Poly (log (./V)) number of elementary gates. For small 
N these networks are very simple. We implemented the measurement of W(a) for N = 4 (two qubits) for 
a variety of initial states. In this case R is a controlled not (CNOT) gate (where the control is in the least 
significant qubit). U is the same CNOT followed by a bit flip in the control. Analogously, V is a sequence of 
controlled phase gates. The complete circuit has at most one Toffoli gate and several two qubit gates. Figure 
3 shows the results of the measurement of the Wigner function for all four computational states of a two qubit 
system. Ideally, W(a) for the state \qo) is nonzero only on the vertical strip at q = 2qo, where it is equal to 
1/2N and on the strip at q = 2qo ± N where it oscillates as (— 1) P /2N. These oscillations correspond to 
interference between the state and its mirror image created by the periodic boundary conditions.^ 

To measure W(a) we used a liquid sample of trichloroethylene dissolved in chloroform. This molecule 
has been used in several three qubit experiments where the proton ( 1 H) and two strongly coupled 13 C nuclei 
{C\ and C2) store the three qubitsP^l We used C\ as our probe particle and the pair H-C2 to store the state 
whose Wigner function we measure. The coupling constants are Jjjc\ = 200.76 Hz, Jhc 2 = 9-12 Hz and 
Jc\Ci = 103.06 Hz while the C\-Ci chemical shift is <5ciC 2 = 908.88 Hz. We determined the value of 
W(a) for each of the independent 16 phase space points. Each of these circuits corresponds to a different 
sequence of r.f. pulses and delays (the number of pulses in each sequence depends on a varying between 
5-17 and taking at most 100ms to execute, -C Xi,T2 of our sample). We used temporal averaging^ to 
obtain, from the part of p that deviates from the identity, the four pseudo-pure initial states whose Wigner 
functions are shown in Figure |3] The experiments were done at room temperature on a standard 500Mhz 
NMR spectrometer (Bruker AM-500 at LANAIS in Buenos Aires and DRX-500 at Los Alamos). We used a 
5mm probe tuned to 13 C and X H frequencies of 125.77 MHz and 500.13 Mhz. The most important sources of 
errors come from the strong coupling and the numerical uncertainty in integrating the spectra. These results 
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illustrate the tomographic measurement of a discrete Wigner function and agree very well with theoretical 
expectation. 

In summary, we presented a general algorithm enabling us to view spectroscopy and tomography as 
dual tasks. From our algorithm follows the construction of circuits characterizing spectral properties of 
general operators and the design of a general tomographer with which we can directly measure the Wigner 
function of the state of a quantum system. The measurement strategy is an example of a very general 
tomographic scheme in terms of which one can interpret previous experiments performed to measure Wigner 
functions. The analogy between this quantum computation and a scattering experiment (that has the dual 
use of providing information either about the state of the scatterer or about the interaction Hamiltonian) is, 
we believe, particularly illuminating and useful. 
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Figure 1: Circuit for measuring Re[Tr(p{7)], for any unitary operator U. 



Figure 2: Circuit for evaluating the spectral density of an operator U. The second register, formed by n\ 
qubits, is prepared in the initial state \E) and is subject to a controlled Fourier transform. The second 
controlled operator acts on the second and third registers mapping the computational states | 2 ® \n)z into 
|t)2®t7*|n)3. 

Figure 3: Measured Wigner functions for the four computational states of a two-qubit system (built with a 
liquid sample of TCE in an NMR spectrometer). Ideally, these Wigner functions should be nonzero only on 
two vertical strips where they take values which are ±1/8. Experimental results show small deviations from 
these values (with a maximum error of 15%. 
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